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THE INFLUENCE OF FREQUENCY UPON THE 
SELF-INDUCTANCE OF COILS. 

By J. G. Coffin. 

Presented by A. G. Webster, January 10, 1906. Keceiyed February 14, 1906. 

It is well known, through the work of Lord Rayleigh,! and Stefan ,2 , 
that currents of high frequency in straight wires distribute themselves 
over the cross-section of the wires very differently from steady currents. 

The values of the resistance of such wires, as well as their self-induc- 
tance, thereby suffer a change. Resistance undergoes a large increase, 
self-inductance a small decrease. The theory for straight wires has been 
hitherto directly applied to wires wound on coils, in want of a better 
theory for such a ease. 

Recently, Max Wien 3 and A. Sommerfeld* have attacked successfully 
the problem for wire coils. 

Max Wien's results, however, are only applicable to small frequencies, 
while those of Sommerfeld, although valid for all frequencies, gave only 
the change in resistance. 

In the following will be found the derivation of a formula by a method 
based on that of Sommerfeld for the change of self-inductance with the 
frequency. 

By means of this theory we shall derive a method of assigning the 
correction to be applied to any coil of a single layer. Its extension to 
coils of more than one layer is reserved for a future communication. 
The theory 5 shows that with increasing frequency the currents tend to 
flow more and more on the inside surface of the coil. The maximum 
change in self-inductance is then the change due to a diminution of the 
mean radius of the coil of one half the diameter of the wire, with which 
the coil is wound. 

1 Phil. Mag., 21, 381 (1886). 
« Wied. Ann., 41, 400 (1890). 

* Max Wien, Ann; d. Phys., 14, 1 (1904). 

* A. Sommerfeld, Ann. d. Phys., 15, 673 (1904). 
8 Max "Wien, loc. cit. 
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For example, consider a coil whose mean radius is 27 cms. wound 
with wire of diameter .06 cms. There is a maximum possible change 
of 3.5 per cent for very high frequency in the self-inductance of a single 

turn, while in a coil of 716 turns 
the total possible change is nearly 
0.3 per cent. This is not a neg- 
ligible amount, so that in standard 
coils and work of precision it is 
necessary to know what the cor- 
rections are for any frequency. 



A' 



tz 



ri 



Derivation op the Formula. 

Consider an infinite circular 
current sheet of thickness d, of 
inside radius r-i and outside r^ . 

Let the symbols S { ) and 
F ( ) denote the scalar and vector 
products respectively. V is the 
Hamiltonian operator. FV ( ) is 
the same as curl ( ) or rot ( ), 
while (SV ( ) is the same as 
div ( ). 

By symmetry, the field inside 
the coil is axial, and a function 
of r alone, at any given time, 
and on account of the coil being 
infinite in length, the field is all 
inside of it. 

In the dielectric, Maxwell's 
equations are, if we neglect the 
displacement current, 

VVH^ and ^ V^ = 0. (1) 



From the second relation we find 
that H is independent of z, and 

from the first relation, that it is independent of x and y ; hence it is 

constant inside the coil, and for 



FiGDBE 1. 



r ^ rj assume H ■= H(,e' 



vat 



(2) 
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where to is 2 w'X frequency, and Hq the maximnm value of H. Outside 
the coil, H must be constant, but as it vanishes at infinity it is zero 
everywhere. Hence, for 

r 5 ra H=: 0. 

In the material of the coil, Maxwell's equations are 

iwq=VVff, ^ = -VVF, a.ndq = <rF, (3) 

from which follow SVB =0 a,ud SVq = (4) 

where q is the current density, 

F is the electric force, 
<r is the specific conductivity. 

Eliminating q and F, by taking the curl of the first equation (3) 

4xFV^= FVrv^= V/SViT- V^^" 

= -V^H by (4) 

and from the second and third of (3) 

^ dt 

therefore 4 jr o- -7- = V^H. (5) 

Now assume ^in the material of the coil to be given by a function Cof r 
alone (on account of symmetry), multiplied by H^e^*, or H= HoUe^* 
where {7=1 when r = rj 

Z7=0 " r = ri. 

The differential equation then which ZZmust satisfy by equation (5) 
becomes in polar co-ordinates 

where k = —i i ir tr <a. 

The solutions of this equation are the Bessel's functions of order zero, 
Jq and K^ and hence 

H= H, {A J, (hr) + BK, (kr)} «««' (7) 

where A and B are arbitrary constants. 
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^— V 7 J)fc(r,-r.) .-i«n-rJ -^» *'"' (^) 



Now, Sommerfeld 8 has shown that the expression 

-r) _ g-<*(r,-r) 
r g«(ra-r,) _ g-»*(»'2-»'i) 

is a very good approximation formula for the case in which It r-i and Ic r^ 
are large, and he finds also that it is valid even when these quantities are 
small, this being an unexpected result. 
Since 

4-g = »^V7 ga-(..-.,_g-..(.-.o -^o^"' (9) 

where in differentiating, on account of the largeness of ^ r only the ex- 
ponential terms were considered variable. 

For the total current /per unit of length we obtain the value 



1= ^f^dr = - -^^-^dr = ^ \Hr, - Hr, 



(10) 



r6 



• T -* J.(i)t 



where i\7' is the number of turns per unit length of the coil. 
The electro-kinetic energy of a system is given by 

T^lzr (11) 

where L is its self-inductance, and / is the total current. In the case 
of alternating currents, the mean square values are meant, or 



Hence we may write 






T 

-J 



Tdt 
Z = ?#=2l±^— (12) 



-:/ 



pdt 





where t is the period. 



« Loc. cjt., p. 678. 
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But the mean square value of an harmonic function is one half the 
square of its amplitude, and in working with imaginary expressions 
the square of the amplitude can always be obtained by multiplying 
any expression by its conjugate. Hence 

T 

p = ^ frdt = I /I (13) 



where /denotes the conjugate to /. 
Therefore 



and 





In general, 

if /A, the magnetic permeability, is assumed equal to 1 everywhere, which 
becomes, in this case, 

CO 

T=^ fH^^Tcrdr 



per unit length of coil 

=:^ jff'rdr + ~ JH^rdr 



Ti 



(16) 



where, in Ti, H=. H^e 

and in T^, H= eq. (8). 

We may write 

i CTdf = ^ frdr^ j'll'dt = ^ Cffffrdr. 



For 71 this becomes ( -j^ j . (17) 

To obtain the integral T2 first put in equation (8) A (1 + i) for tk. 
Then, calling h (r^ — r)=y, 
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and Ji {r.^ ~ n) = y^, 

we obtain for H&uA its conjugate H, the expressions ; 

/r, piy pV — fi—^v p—y 

V J. g»^i gVi _ g-tyi g-vi 

and (18) 



Jn e-^y eV - e'y e-y ^ ^_^^ 






A conjugate expression is obtained by putting — i for i everywhere in 
an equation. 

Their product becomes 

or 

^jy_ jys^ cosh2y-cos22/ 



r cosh 2yi — cos 2 «/x 
Multiplying this expression by r dr and integrating from r-i to r^ gives 



/ 



cosh r dr ■=■ sinh r 

putting a; = 2 ^ (ra — >-i) = 2 A c? 

^0^ ''i <^ sinh a; — sin a; 

" 8 a; cosh a; — cos x 
"We may now write 

where x — ld V2 ir w o- = (4 ir c? Vct^) V** where w = 2 tt w 

Now, 4 TT^iV^^j-i^ is the self-inductance per unit length of a coil with a 

mean radius r-i, and is therefore the self-inductance for infinite frequency, 

according to the remarks of Wien. This is also shown by equation (19), 

sinh X 
as the second term in parenthesis is zero when a; = co for — -r— = 1, 

and sinh x and cosh x increase without limit as x increases without limit, 
as may be seen by their expansions 



3? X^ 

sinh a; = a; + jr-7 + ^n+. . .;sina; = 
3 ! ! 


X* x- 


cosh a;=14- — + — + . . .;cosa; = 


2!^4! 


juation (19) may then be written 




^» — ^00 _2 c^ /3 sinh a; - 


sin kN 2 <i 
coax) ~ 3 j-j 
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(20) 



(21) 

By this equation we see that the fractional change of self-inductance in 
any coil may be calculated if the properties of the function 

^ ^ 3 sinh a: - sin a; 
X cosh X — cos X 

be known for all values of x from zero to infinity. 

A formula for small frequencies may be derived by expanding the 
hyperbolic and circular functions and retaining only low powers of the 
variable x. Q then becomes 

la;* a;' 

3I + 7I + Yv + ■ ■ * 
9 = 3 



1 a;* x^ 



10! 
which, retaining only terms in x^ is 

1 3 a;* 6 a;* 

^ 2"^TT ^ "^TT _. ^_ 14£* _ , _ 8^* _ -, _ ^ 
'i^i: '1^^ ^' 7!~ 7!" 630' 

2 ■*" a;* "•■ 6 ! 

This equation is applicable for small frequencies only, and (21) may be 
written, since 64:7^ = 630 very approximately. 

It agrees as to its form and dimensions with that of Max Wien.7 When 
<o = 0, it reduces to 2/3 d/ri. As a verification of the term 2/3 d/r^, 
the approximation formula for Bessel's formula with small argument 

H=H,diri:ZJLe^*^ (24) 

y r r^ — r^ ^ ■' 



' Loc. cit., p. 16. 8 Soramerfeld, loc. cit, p. 679. 
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was taken and a calculation similar to the one above for the general case 
made. This completely verified this term, and proved that the expres- 
sion used above equation (21) is actually valid for small arguments as well 
as large. This last expression 2/3 d/r^ multiplied by Zoo = iiir^ N'^ r^ is 

Xo-Zoo =|^''iV^'ritf (25) 

o 

and is therefore the additional self-inductance due to the field in the 
wires themselves for steady values of the current, or, in other words, 
this is the maximum possible change in self-inductance. Calculations by 
means of (25) show this expression to give, for coils of finite length and 
wires with round instead of square cross-sections, the values : 





TABLE 


I. 




PARISON OP 


THE Theory with the 


Tkue Change in Self-I 


NDUCTA 


Coil Length. 


True Change. 


Value of 2/3 «f/ri. 


Ratio. 


46.0 cms. 


.000380 Henry 


.000280 Henry 


1.35 


30.5 « 


.000200 " 


.000160 « 


1.25 


15.0 « 


.000067 « 


.000057 « 


1.17 


13.0 « 


.000051 « 


.000044 « 


1.16 



The column headed " True Change" was computed by taking the differ- 
ence between the self-inductance of a coil of mean radius a, computed by 
means of an exact formula, and that of the same coil, using as a mean 
radius a — p, where p is the radius of the wire. These results, consider- 
ing the assumptions made in the theory, are in good agreement with it. 
To make the theory fit to actual facts the second term in the parenthesis 
in equation (19) should be multiplied by a constant, the average value 
of which is approximately 1.25, deduced from the above table. We 
should expect some multiplier to be necessary because a current sheet is 
the equivalent of square wires without insulation, while, in the coil for 
which the calculations were made the wires are round and have insulation. 
We may then write 



^"~^°° = 1.25 



2 d f Z sinh a; — sin a; \ 

3 rx \x cosh X — cos xj' 



In the proposed use of this theory, however, it is not necessary to em- 
ploy this constant whose value at best is not satisfactorily determinate. 

As cosh X and sinh x increase without limit as x increases, equation 
(26) may be written, for large values of <«, 
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L^ — L^ , „. 2 rf 3 b d 



Znx ^nd'\/2-7ru>cr 



797 
(27) 



This formula agrees to a constant factor pres, with that of Lord Ray- 
leigh,^ for large values of n for straight wires ; i. e. : 

Zoo Y 2we dy'iiTriav 

TABLE II. 

Table as Valttbs tor the Curves Q, Qj, and Q^, fob Values of x from 

Zero to Infinity. 



X. 


Gosh X. 


Sinh X. 


Cosx. 


Sinx. 


Q- 


ft- 


ft- 


0.0 


1.0000 


0.0000 


1.0000 


0.0000 


1.0000 


1.0000 


00 


0.5 


1.1276 


0.5211 


0.8776 


0.4789 


0.9999 


0.9999 


6.0000 


1.0 


1.5431 


1.1752 


5402 


0.8415 


0.9984 


0.9984 


3.0000 


1.5 


2.3623 


2.1291 


0.0706 


0.9975 


0.9920 


0.9918 


2.0000 


2.0 


3.7622 


8.6269 


-0.4163 


0.9092 


0.9756 


0.973 


1.5000 


3.0 


10.0677 


10.0178 


-0.9900 


0.1409 


0.8932 


0.871 


1.0000 


4.0 


27.308 


27.262 


-0.6634 


-0.7570 


0.7516 


0.694 


0.7500 


5.0 


74.205 


74.197 


+0.2840 


-0.9588 


0.6100 


Negative 


0.6000 


6.0 


202 approx. 


202 approx. 


0.9603 


-0.2790 


0.5032 




0.5000 


7.0 


573 approx. 


573 approx. 


0.7537 


+0.6574 


0.4287 




0.4289 


8.0 

9.0 

10.0 

00 




Very large 

and the same 

as sinh x. 


- 


Very large 

and the same 

as cosh X. 






0.3750 
0.33-33 
0.3000 
0.0000 




0.3750 
0.3333 
0.3000 
0.0000 


CO 


00 


00 


Indeterm. 


Indeterm. 





—00 






Computations were carefully made for the values of the term 

- 3 sinh X — sin a; 

X cosh X — cos X 



(28) 



' Loc. cit. 
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for different values of x, as well as of ^i and ^2; where 

^^ = (l - ^) ^"'l (29) 

Q, = \. (30) 

Table II contains all of the data. 

It is clearly seen that as long as a; < 2 that formula ^i may be 
used, while if a; > 7, formula Q^ is valid. In between these points, the 
calculated values will serve. 

A plot of these three functions on the next page shows still more 
clearly the relations of the three formulae for different values of x. 

Manner of Using the Curve. 

This curve is to be used in the following manner. We assume that 
in any coil of a single layer, the falling off of the self-inductance takes 
place according to equation (22), represented by the full-line curve, and 
the maximum ordinate unity is taken to be the maximum change possi- 
ble of the self-inductance, which is otherwise calculable. For any given 
dimensions of the coil, the values of x correspond to definite frequencies, 
and either a table or a plot of corrections may be derived from the above 
curve. 

For small values of x, and therefore for small values of the frequency, 
the correction varies proportionately to d^, while for very large frequen- 
cies it is proportionate to \/d by equations (29) and (30). 

As the first case is the more usual, it is important, in order to avoid 
the necessity of any correction, to wind the coil with wires of small 
diameter. The size necessary to avoid a correction under a given fre- 
quency may easily be calculated. As an application of this theory, the 
following example will be carried out: 

Consider a coil of 27 cms. radius = r-^ wound with a single layer of 
wire whose diameter d is .063 cms. Let the conductivity of the wire be 
taken as 5.9 X 10"^. Then the relation between x and w is w = 2630 a:^ 
from (19) where n = frequency. In the following table, column 2, are 
the values of x corresponding to the frequencies n in column 1. In 
column 3 are the values of the ordinates corresponding to these values of 
X taken from the Plot (I). 
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PLOT I. 




„ -. 3 smh a; — sin a; , „ ,. 

Plot of the functions Q = r luil line. 

X cosh X —QOi X 



9 10 U 12 13 11 15 
Qi = (1 - x</630) (lotted line. 
Q2 ~ 3/^ (lotted line. 
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TABLE III. 



n 


X 


Ordinates. 








1.000 


500 


0.436 


1.000 


1000 


0.616 


0.999 


1000 


0.755 


0.998 


3000 


1.068 


0.996 


4000 


1.230 


0.995 


5000 


1.380 


0.993 


6000 


1.518 


0.990 


10000 


1.950 


0.977 


15000 


2.380 


0.950 


20000 


2.7G0 


6.916 


30000 


3.880 


0.846 


40000 


3.900 


8.765 


50000 


4.350 


8.700 


60000 


4.780 


0.638 


70000 


5.150 


0.590 


8C000 


5.510 


0.550 


90000 


5.850 


0.516 


100000 


6.170 


0.483 


120000 


6.760 


0.463 


150000 


7.550 


0.397 


170009 


8.050 


0.373 


200000 


8.720 


0.344 



The plot of the corresponding value of n as abscissae and values of 
column 3 as ordinates gives a curve of corrections for the above coil. 
For example, at a frequency of 50,000 there is 0.7 of the total possible 
change (equal to Zg 2/3 d/r-i, X 1.25 approximately) still to go. Or, the 
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PLOT II. 




S 



in "^ 
o a> 00 b- <o >o -^ ^c53 

1-4 O O O O O O 00-~r- 

Curve showing the fractional part of the maximum decrease for different fre- 
quencies in a coil, the wire of which has a diameter of .063 cm. 
Ordlnates — fraction of total change left unchanged. 

. , , frequencies 
Abscissae = —j^gQ 

For example, at 30,000 /-»^ per second the correction is 0.157 X 0.0000618 Henrys 
in 0.0362 Henrys, or about 1 part in 3600. 

VOL. XLI. — 51 
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correction is 0.3 of the maximum change. The value of the correction 
may be best obtained by reading down from the ordinate unity. 

These computations also show that equation (29) 

X* 

may be employed up to frequencies of 10,000 for this particular coil, and 
that equation (30) 

3 
^^ = ^ 

may be used for frequencies beyond 120,000 per second. As the max- 
imum correction is 0.2 per cent of the total, if frequencies under 4000 per 
second are used, the error will not be over 0.2 per cent X 0.005 = 0.001 
per cent, which is practically negligible. 

If, however, other things being equal, the diameter of the wire were 
four times larger than it is, the same corrections would correspond to 
frequencies 16 times less, or an error of 0.001 per cent is now obtained 
at a frequency of only 250. 

This suggests, if large currents must be carried, that to wind the coil 
with flat wires would diminish the effect of the frequency upon both the 
changes in self-inductance and resistance. 

The Plots I and II, just derived, show how the self-inductance of a 
circular current sheet falls off with the frequency. In other words, the 
coil is assumed to be wound with wires of rectangular cross-section. 
This may cause a doubt as to the applicability of the curve in question. 

In the first place, the method of its use does away with the necessity 
of knowing the constant multiplier of the function Q, as explained above ; 
and in the second place, the following two reasons go to show the entire 
correctness of the deductions. 

Sommerfeld finds that his results agree to a constant pres with those 
of Max Wien, which were deduced for circular wires, and both formulae 
agree remarkably well with the experimental data at hand ; this for the 
increase of resistance with the frequency. We should, therefore, fairly 
expect the results for the decrease in self-inductance to also be correct to 
within a constant factor. As there are no data available for the decrease 
in self-inductance, we have not been able to verify this conclusion. To 
further assure ourselves, however, we have worked out the problem of 
the change of self-inductance with the frequency in an infinite straight 
circular-ring conductor, using the approximation formulae for the Bessel's 



COFFIN. — THE SELF-INDUCTANCE OP COILS. 803 

functions as above given. It was found that the variable part of the 
expression for self-inductance agreed identically with that derived for a 
circular current sheet. This shows that the bending of straight wires 
into circles, and forming a coil of a single layer therewith, does not affect 
the wai/ in which the self-inductance varies. 

We believe, then, that Plot I gives very approximately the law 
according to which the self-inductance falls off for straight ring-shaped 
conductors, and probably for any simply shaped square or polygonal 
hollow straight conductors, as well as for circular coils of a single layer, 
wound with wires of circular, square, or polygonal cross-sections. 

Dept. Physics, 
Claek University. 



